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The recently discovered three-dimensional hyperhoneycomb iridate, has raised hopes for the real¬ 

ization of dominant Kitaev interaction between spin-orbit entangled local moments due to its near-ideal lattice 
structure. If true, this material may lie close to the sought-after quantum spin liquid phase in three dimensions. 
Utilizing ab-initio electronic structure calculations, we first show that the spin-orbit entangled basis, jeff =1/2, 
correctly captures the low energy electronic structure. The effective spin model derived in the strong coupling 
limit supplemented by the ab-initio results is shown to be dominated by the Kitaev interaction. We demonstrated 
that the possible range of parameters is consistent with a non-coplanar spiral magnetic order found in a recent 
experiment. All of these analyses suggest that may be the closest among known materials to the 

Kitaev spin liquid regime. 


Introduction - Kitaev’s exact solution of a quantum spin- 
liquid on a spin-1/2 honeycomb model has spurred consider¬ 
able interest in the search for a material realization[l, 2]. Of 
particular focus is the family of quasi-two-dimensional (2D) 
honeycomb iridate materials Cf-A 2 lr 03 {A = Na, Li, here¬ 
after aAlO), where iridium (Ir) ions form decoupled layers of 
honeycomb lattices [3, 4] and have been argued to host spin- 
orbital entangled jeff = 1/2 degrees of freedom[5-7]. Due 
to the interplay of strong atomic spin-orbit coupling (SOC) 
and correlation effects, these jeff = 1/2 moments in the ideal 
aAlO structure interact in the highly anisotropic manner de¬ 
scribed by the Kitaev model[8]. In addition to these Kitaev- 
type exchanges, the symmetries of the ideal structure also per¬ 
mit additional exchanges that generate a plethora of interest¬ 
ing phases of matter[9]. In reality, however, these materials 
possess sizeable monoclinic distortions that deform the octa¬ 
hedral oxygen cages surrounding Ir ions [4, 10]. These distor¬ 
tions lower the symmetry of the system and therefore com¬ 
plicate the description of these materials. Thus far, a consen¬ 
sus on the minimal model required to describe this family of 
2D honeycomb iridates has yet been reached; a distortion-free 
analog of these honeycomb iridates may offer a more direct 
path towards the realization of Kitaev physics. 

The timely discovery and synthesis[12, 13] of the hyper¬ 
honeycomb / 5 -Li 2 lr 03 (hereafter /3LIO) may present such an 
exciting opportunity. Much like its 2D counterpart, the Ki¬ 
taev model on the ideal, 3D hyperhoneycomb lattice supports 
an exact spin liquid ground state[ 14-17]. In addition, the 
distortion-free, classical pseudospin-1/2 model on the hyper¬ 
honeycomb lattice also supports a myriad of complex mag¬ 
netic phases[ll]. Moreover, interesting topological phases 
have been predicted on this lattice[18]. These previous re¬ 
sults illustrate the possibilities that may be realizable in /3LIO; 
however, they rely on the use of the jeff = 1/2 degrees of 
freedom in the low-energy description of /3LIO, which has 
not been justified microscopically. Furthermore, whether the 
near-ideal structure of /3LIO can give rise to a simple mini¬ 
mal pseudospin model dominated by the Kitaev exchange has 
so far not been validated. Also, with the recent experimental 



FIG. 1. (Color online) Phase diagram of the J-K-V model repro¬ 
duced from Ref. [ i 1 ], overlaid with density distributions of exchange 
interaction parameters estimated from ab-initio results for /3LIO; see 
main text for details. Six shaded areas with different colors represent 
relevant Uqe values ranging from 1.5 to 3.0 eV in increments of 0.3 
eV. Vertical and horizontal markers on the shaded areas depict the 
variation of Uqe and Jh, respectively, as shown in the inset. Note 
that, when Uqe > 2.4 eV, corresponding phases lie within the green 
SPa- spiral phase area. SPa- phase is consistent with the magnetic 
order observed in the experimental work of Ref. [12]. For detailed 
discussion of the other phases, see Ref. [11]. 


observation of a magnetic spiral order in /3LIO, any minimal 
model and its accompanying parameters must also be capable 
of predicting the observed order: this provides a stringent test 
of feasibility for any model describing /3LIO. 

In this letter, we tackle these issues by combining results 
of our ab-initio electronic structure calculations and a strong¬ 
coupling theory to arrive at a jeff = 1/2 model to describe 
/3LIO. From our ab-initio band structure results, we find that 
the low-energy states can be described in terms of localized 
Jeff = 1/2 states because of the large atomic SOC present in 
Ir. In fact, the magnitude of SOC in the paramagnetic state is 
enhanced by the electron interactions in Ir d orbitals, which is 
consistent with recent observations in several Ad and 5d transi- 
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FIG. 2. (Color online) (a) Network of IrOe octahedra in the hyper¬ 
honeycomb lattice. The two distinct nearest-neighbor (NN) bonds, 
X and Z, are depicted as solid green and red lines, respectively. X 
and X’ bonds are symmetry equivalent, whereas Z bonds are dis¬ 
tinct. Emerging from the octahedra, the red and blue arrows point 
in the direction of the trigonal distortions for red and blue IrOe oc¬ 
tahedra respectively. The trigonal distortion, which consists of the 
compression and rotations of the opposing oxygen trangles, is illus¬ 
trated in the figure and in the inset, (b) shows the band structure with 
the presence of spin-orbit coupling (SOC). Solid red and dashed grey 
curves are the band structure of the experimental and ideal structures, 
respectively. 


tion metal compounds[ 19-21]. To go beyond the limitation of 
ab-initio calculations in treating electron interactions, we em¬ 
ploy the strong-coupling expansion recently proposed in [9] to 
arrive at a minimal pseudospin-1/2 model. Due to the near¬ 
ideal structure of /3LIO, we discover that the resulting pseu¬ 
dospin model is near-isotropic while both distortion-induced 
and further neighbor interactions are small. Remarkably, the 
estimated exchange interactions places the model near the fer¬ 
romagnetic Kitaev limit and within a region where the classi¬ 
cal ground state agrees well with the experimentally-observed 
spiral phase[12]. Our results are summarized in Fig. 1 and 
elaborated in the rest of this work. 

Structure and ab-initio calculations - /3LIO is a mem¬ 
ber of the generic three-dimensional (3D) harmonic honey¬ 
comb iridate series [22], which are structural variants of 2D- 


honeycomb iridates aAlO. The hyperhoneycomb lattice is 
composed of a tri-coordinated network of edge-shared IrOe 
octahedra as shown in Fig. 2(a). There are two types of 
nearest-neighbor (NN) bonds in this network: we denote 
these bonds as X and Z. Despite being symmetry-inequivalent, 
these two bonds are almost identical owing to their simi¬ 
lar local crystal structures as revealed by recent structural 
analysis [12, 13]. 

The crystal structure refinement also revealed nearly ideal 
IrOe octahedra compared to those from aLIO: standard de¬ 
viations of the Ir-0 bond length and O-Ir-0 bond angles in 
^LIO (0.002 A and 2.11°) are much smaller than those in 
OfLIO (0.050 A and 4.56°)[10]. Since finite standard devia¬ 
tions in O-Ir-0 bond angles is a result of trigonal distortion in 
the IrOe octahedra, the small value present in /3LIO indicates 
that trigonal distortions are indeed small in this compound (the 
directions of trigonal distortion are shown as colored arrows 
in the figure). The nearly ideal IrOe octahedra in /3LIO sug¬ 
gest that local crystal fields are principally cubic in symmetry, 
therefore the spin-orbital entangled jeff = 1/2 states would 
be a good basis to construct a low-energy description of this 
material in the presence of strong SOC. 

To validate the use of jeff = 1/2 states in the low-energy 
description of /3LIO, we turn to ab-initio electronic structure 
calculations [23]. The band dispersions of the ideal and exper¬ 
imental structures with SOC can be seen in Fig. 2(b). The 
dispersions from the experimental structure (solid red curves) 
and those from ideal one (dashed grey curves) share similar 
overall shape, especially near the chemical potential. The 
separation between the upper eight bands and the lower six¬ 
teen bands (including Kramers degeneracies) can be clearly 
seen in the figure, suggesting the formation of jeff = 1/2 and 
Jeff = 3/2 bands [24]. 

In Fig. 3(a), we show the jeff -projected band dispersions 
and density of state (PDOS) in the presence of SOC based on 
the experimental structure. The Projection is done by taking 
inner products between the atomic jeff states and the Bloch 
state represented in terms of the local pseudo-atomic orbital 
basis. Weights of the jeff =1/2 and 3/2 components within each 
Bloch state are depicted as the size of red and blue circles, re¬ 
spectively, in the band plots. The large jeff = 1/2 PDOS 
weight in the upper eight bands—the closest bands to the 
Fermi level—indicates the development of jeff =1/2 bands 
and confirms that the basis states relevant to the low-energy 
description of /3LIO possess mostly jeff =1/2 character. The 
effect of electron correlations inherent to Ir t 2 g orbital further 
enhances the jeff = 1/2 character as shown in Fig. 3(b), where 
effective on-site Coulomb interaction = U — Jh in¬ 
cluded within the DFT-^f/ formalism {Jh is Hund’s coupling; 
for details see Sec. A in Supplementary Materials). [25] As 
t/eff = 3.0 eV is added, the separation between the jeff = 1/2 
and Jeff = 3/2 states becomes enlarged. This SOC enhance¬ 
ment is also manifested in the increased magnitude of the ef¬ 
fective t 2 g SOC At 2 g as shown in Table I, which is obtained 
from the on-site matrix elements in the Wannier orbital calcu¬ 
lations [26]. Such behavior has also been reported in other 4d 
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FIG. 3. (Color online) (a,b) Band structure and density of states (DOS) projected onto the jeff states in the presence of SOC (a) without and 
(b) with the on-site Coulomb interaction f/gff = 3.0 eV. (c) shows the schematic shape of the jeff=l/2-like Wannier orbital constructed from the 
ieff=l/2-dominated bands near the Fermi level. Dashed and solid circles depict the Wannier orbitals from calculations without and with finite 
f/eff, respectively. Weights of the central jeff=l/2 and nearest-neighboring jeff=3/2 tail in the orbital are shown in (d) as a function of f/eff. 


t/eff 

(eV) 

0.0 

1.5 

3.0 

At2g 


0.401 

0.482 

0.516 

tl 

z 

+0.085 

+0.077 

+0.064 


X 

+0.083 

+0.074 

+0.058 

|t2| 

z 

0.238 

0.255 

0.270 


X 

0.260 

0.276 

0.289 

ts 

z 

-0.162 

-0.119 

-0.060 


X 

-0.153 

-0.110 

-0.055 


TABLE I. Magnitude of SOC within the Ir t 2 g states and t 2 g hopping 
terms from Wannier orbital calculations in the presence of t/eff • We 
adopt the coordinate system such that t 2 is negative for both Z and X 
bonds. By symmetry, t 2 is positive for the X’ bonds. 


and 5d orbital systems [27][19-21]. 

The effects caused by electron correlations in the low- 
energy Jeff = 1 /2-dominated states deserve further comment. 
Fig. 3(c) shows the schematic shape of Wannier orbitals con¬ 
structed from the jeff = 1/2 energy window, which can be 
considered as the local orbitals that span the low-energy sub¬ 
space. Owing to the nearly ideal IrOe octahedra (as sup¬ 
ported by the small amount of trigonal distortion of less than 
100 meV), the Wannier orbitals consist of pure jeff = 1/2 
character on the center Ir site, while it has jeff = 3/2 tails on 
the NN sites. Similar features have been reported in cfNIO and 
CfLIO[28, 29], which mirrors the remnant molecular orbital 
character originating from the t 2 g hopping [30]. As f/eff is in¬ 
cluded and At 2 g is enhanced, the jeff =1/2 character becomes 
more dominant while jeff = 3/2 components on the NN sites 
decreases as shown in Fig. 3(c) and (d). The jeff = 1/2-like 
Wannier orbital is more localized accordingly, which makes 
the low-energy description of /3LIO in terms of the localized 
Jeff = 1/2 states more feasible in the strong coupling limit. 

t 2 g Wannier orbital hopping amplitudes - For a detailed un¬ 
derstanding of how the near-ideal structure of /JLIO is man¬ 


ifested in the electronic band structure, we calculated the Ir 
t 2 g hopping amplitudes from the Wannier orbitals in the ex¬ 
perimental structure. Table I shows the magnitude of the three 
largest hopping terms—ti, ^ 2 , and ts —as the value of (7eff 
changes ((7eff = 0.0 eV, 1.5 eV, and 3.0 eV, and SOC is in¬ 
cluded in the calculation); see Sec. B and Fig. Sl(a) in Sup¬ 
plementary Materials and Ref. [31] for illustration of these 
hopping processes. Since the Ir-Ir bond lengths and Ir-0- 
Ir bond angles are similar on the two inequivalent bonds of 
/JLIO (X and Z bonds), the values of their respective hopping 
amplitudes are expected to be similar. Indeed, by compar¬ 
ing the hopping amplitudes between the two inequivalent NN 
bonds, we observe small anisotropies between the X and Z 
bonds (< 10%), which reflects the close-to-ideal structure of 

puo. 

The evolution of the NN hopping amplitudes as we include 
on-site Coulomb interactions can be seen in Table I. As Ceff 
increases, |t 2 | increases while ti and ts decrease. Such be¬ 
havior is understood in terms of the enhanced hybridization 
between the Ir t 2 g and oxygen p states in the presence of (7eff. 
Inclusion of Ceff pushes the jeff=3/2 states down energetically 
so that they become closer to the oxygen p states. This leads 
to increased hybridization between the Ir t 2 g and oxygen p 
states, which yields the enhancement of oxygen-mediated t 2 
(and the reduction of and ts). 

Strong-coupling minimal model and experimental spiral 
phase - Having validated the use of the jeff = 1/2 basis 
and the similarity of hopping amplitudes between inequiva¬ 
lent bonds, we can now construct an effective model to de¬ 
scribe the low-energy properties of /JLIO in the large-(7 limit. 
Following the derivation in Ref. [9], we start with localized 
Jeff = 1/2 states then perform a strong-coupling expansion us¬ 
ing NN t 2 g hopping amplitudes [32]. In the presence of Hund’s 
coupling Jh, we arrive at a NN, jeff =1/2 model with highly 
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anisotropic pseudospin exchanges 

H=Y^ rSi-Sj + K°‘S^S‘}+ T°‘{S^S]+ S'lS'^), 

{ij)ea{M) 

where Si is the jeff = 1/2 pseudospin on site i, a labels the 
NN {ij) bond by its Kitaev component, and [3 and 7 denote 
the two non-Kitaev components of the (ij)-bond. The ex¬ 
changes J, K, and V are functions of the hopping amplitudes 
ti-ts, strength of Hund’s coupling Jh, SOC A, and the on-site 
Coulomb interaction U : the relation between these quantities 
are given in Sec. D in Supplementary Materials. 

To establish the region in the parameter space that best 
models ^LIO, the following statistical analysis was employed. 
First, the hopping amplitudes and SOC values in Table I 
were interpolated as a function of f/eff. Next, f/eff and Jr 
were treated as independent variables with choices of ranges 
1.5 < f/eff < 3.0 eV and \t 2 \ < Ju < Xt 2 g (note that ^2 is 
the largest hopping term)[33], due to the difficulty in deter¬ 
mining specific values of f/eff and Jr. In order to present the 
phase evolution as a function of Ceff ^ we chose six f/eff inver- 
vals centered at 1.5, 1.8, 2.1, 2.4, 2.7, and 3.0 eV with ranges 
Af/eff = iO-1 eV. From these parameters, six possible ranges 
of the exchange parameters were estimated as shown in the 
inset of Fig. 1[34]. We found that the mean anisotropies be¬ 
tween X and Z bonds in J, K, and F are 3%, 15%, and < 1% 
respectively relative to the largest energy scale, which is the 
Kitaev exchange. As a first approximation, we treated all ex¬ 
changes as isotropic between the X and Z bonds, which yields 
the NN Hamiltonian studied in the classical limit in Ref. [11]. 
Lastly, we overlaied the density distributions of the exchange 
parameters on top of the relevant portion of the classical phase 
diagram reproduced from Ref. [11], thereby yielding Fig. 1. 

The phase diagram is the quarter of a polar plot near the 
ferromagnetic Kitaev limit: The angular coordinate (j) shows 
the ratio between J and K — tan((/)) = J/K. Meanwhile, 
the radial coordinate depicts the strength of F — 0 G ( 7 r/ 2 , tt) , 
tan 6 > = y/ J‘^ K‘^/T. The bottom boundary of the half-ring 

(6> = 7 r/ 2 , 0 = 37 r/ 2 ) is the Heisenberg-Kitaev limit (F = 0) 
and the origin (0 = tt ) is the pure F limit. As seen in the 
figure, the NN exchanges in f3LlO are likely dominated by a 
large, ferromagnetic Kitaev exchange, perturbed by small \J\ 
and |r|. Note that, higher t/eff and Jr prefer smaller 0 (larger 
K) and smaller 0, respectively. 

We find that, when f/eff >2.4 eV, the correspondins phases 
lie within area of the spiral phase SP^ . Remarkably, this non- 
coplanar spiral magnetic phase possesses the same symme¬ 
tries as the experimentally determined magnetic order.[l 1 , 12 ] 
In other words, using the ab-initio hopping and SOC param¬ 
eters, the resulting exchange parameters in the isotropic J- 
K-T pseudospin model results in a classical ground state that 
agrees with the experimental magnetic order. 

Discussion and Conclusion - Although we have shown that 
the Jeff = 1/2 states form a valid basis as a consequence of the 
small amount of distortions present, the difference between 
the dispersions of the ideal and experimental structures away 


from the Fermi level are due to these distortions and the result¬ 
ing bond anisotropies. However, in the context of the effective 
pseudospin model, we have shown that these non-idealities 
are negligible for the J and F exchanges. The Kitaev ex¬ 
change, on the other hand, is more anisotropic between the 
X and Z bonds, but we speculate that the SP^- spiral phase 
will remain robust under this anisotropy; we leave the investi¬ 
gation on the effects of bond anisotropy for future work. Nev¬ 
ertheless, it remains true that the NN exchanges on both X 
and Z bonds are dominated by large ferromagnetic Kitaev ex¬ 
changes and that the Kitaev spin liquid is robust against bond- 
anisotropies.[l] 

In addition to distortions and bond anisotropies, an accu¬ 
rate description of the electronic structure also requires hop¬ 
ping amplitudes beyond the NN level (see Sec. C in Supple¬ 
mentary Materials for details). These terms would generate 
further neighbor exchange interactions in the strong coupling 
theory. However, these exchanges are no more than 10% of 
those at the NN level. Since we expect that such small further- 
neighbor interactions do not change our conclusions, we fo¬ 
cused on the NN exchange interactions in out manuscript. 

To enhance the Kitaev exchange relative to other interac¬ 
tions and to approach the spin-liquid regime of the Kitaev 
model, strengthening the oxygen-mediated-type hopping (^ 2 ) 
is a viable option. Increasing the on-site Coulomb interaction 
can further localize the t 2 g orbitals, which reduces the ampli¬ 
tudes for direct hopping channels like ti and ts while oxygen- 
mediated hopping channels like ^2 are comparatively less af¬ 
fected. In addition, increasing U has the effect of driving the 
system deeper into the Mott insulating regime and reducing 
the strength of further neighbor interactions. Therefore, a 4J 
variant of the p\AO may offer the right ingredients to enhance 
the Kitaev exchange. 

Indeed, the isoelectronic, 2D honeycomb Qf-Li 2 Rh 03 has 
been synthesized and argued to be a relativistic Mott insu¬ 
lator driven by electronic correlations and SOC. [35] Further¬ 
more, this material does not magnetically order down to 0.5 K, 
which is an indication of magnetic frustration. [35] We spec¬ 
ulate that the hypothetical 3D polymorph—hyperhoneycomb 
/ 5 -Li 2 Rh 03 —may be a less distorted version of a-Li 2 Rh 03 
that has all the right properties to further approach the Kitaev 
region, in analogy to /3LIO. 
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FIG. SI. (Color online) (a) Ir t 2 g Wannier orbitals for up spin component located at an Ir site, obtained in a GGA+SOC calculation. Other 21 
t 2 g Wannier orbitals can be obtained by translating them to other Ir sites in the unit cell and applying time-reversal operation. Only nearest- 
neighboring Ir atoms are shown, (b) One Ir jeff=l/2 Wannier orbital obtained in a PBE-hSO-i-f/ calculation with f/eff = 3 eV. Other eight 
ieff=1/2 Wannier orbitals can be obtained by translations and time-reversal operation. Note that, isosurface value for drawing the t 2 g Wannier 
orbital in (a) is twice larger than the one used in (b). 


Details on ab-initio electronic structure calculations 

For the electronic structure calculations with SOC and 
on-site Coulomb interaction, OPENMX code[36], which is 
based on the linear-combination-of-pseudo-atomic-orbital ba¬ 
sis formalism, was used. A non-collinear DFT scheme and 
a fully relativistic j-dependent pseudopotential were used 
to treat SOC, and Perdew-Burke-Ernzerhof parametrization 
of the generalized gradient approximation (GGA) was cho¬ 
sen for the exchange-correlation functional[37], which was 
compared and found to be almost identical with the results 
with the Perdew and Zunger local density approximation 
functional[38]. 400 Ry of energy cutoff was used for the real- 
space sampling, and 9x9x9 k-gvid was adopted for the 
primitive unit cell. Electron interactions are treated as on¬ 
site Coulomb interactions via a simplified hDA-\-U formal¬ 
ism implemented in OPENMX code [39], and up to 3.0 eV of 
f/eff = U — Jn parameter (Jh is Hund’s coupling) was used 
for Ir d orbital in our GGA-i-SOC+t/ calculations. Maximally- 
localized Wannier orbital method[40], which is implemented 
in OPENMX code[41], were used to obtain the tight-binding 
Hamiltonian for Ir t 2 g atoms. 

Supplementary Material B: 

Ir t 2 g and jeff=1/2 Wannier orbitals 

In order to obtain the hopping integrals between the Ir t 2 g 
states, the Wannier orbitals for the Ir t 2 g bands were calculated 
in the presence of SOC. Eig. SI shows the results, where the 
three / 2 g Wannier orbitals at an Ir site with spin up component 
are shown. Other 21 orbitals at the four Ir sites in the unit cell 


are obtained by translating them and applying time-reversal 
operation icr^/C, where ay and JC are the Pauli matrix acting 
on the spin sector and the complex conjugation, respectively. 
The oxygen p-orbital components hybridized into the Ir t 2 g 
bands are manifested as the oxygen p-orbital tails shown in 
the figure. Size of oxygen hybridization is slightly enhanced 
as the value of f/eff is increased, which contribute to the ^2 
hopping term dominated by the oxygen-mediated channel. 

The values of the effective SOC strength mentioned 
in the main text were estimated from the on-site energies be¬ 
tween the t 2 g Wannier orbitals. The on-site energy matrix can 
be approximately expressed as i^on ~ ^t 2 g^t 2 g * ^ where 

At is the minor trigonal crystal field terms. The value of Xt^g 
for each Ceff value was obtained by taking the average of the 
matrix elements corresponding to the SOC term, with their 
standard deviation smaller than 10% of their average. 

Pig. Sl(b) shows the jeff=1/2-like Wannier orbital, ob¬ 
tained from the low-energy window dominated by the jfeff=l/2 
character as illustrated in Pig. 3(c) in the main text, in a 
GGA-^SOC+f/ calculation with t/eff = 3 eV. Like the t 2 g or¬ 
bitals, other seven jieff=1/2-like Wannier orbitals in the unit 
cell can be obtained by translation and time reversal opera¬ 
tions. Decomposing the Wannier orbital in terms of the lo¬ 
cal Ir t 2 g basis reveals the dominant jfeff=l/2 character at the 
center with the jeff=3/2-dominated tails on the three next- 
neighboring Ir sites, as schematically shown in Pig. 3(c) in the 
main text. The tail components are gradually reduced as the 
f/eff value is increased, so that the Wannier orbital becomes 
more localized in the presence of higher ^4^. 

Supplementary Material C: 
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FIG. S2. (Color online) (a) Ir sites projected onto xy-plane defined in terms of local cubic axes, ai, a 2 and as are the Bravais lattice vector 
for the primitive cell, where four sublattices within the primitive cell are labeled as Irl to Ir4. Centers for rotations and inversion are 

depicted as red square and black dots on the Z and X (X’) bonds, respectively. Note that d becomes same with the Ir-0 bond length in the 
absence of trigonal distortion. Inset shows the three major hopping channels between NN Ir sites. (b,c) All NNN and (d) third NN neighbors 
for an Ir site, Ir4, are depicted as colored arrows. Note that different colors (blue and green) are used to distinguish Ir sites and bonds belonging 
to one zigzag chain to another, (a) shows NNN hopping paths that can be reached through intermediate NN Ir sites. Intra- and inter-chain 
bonds are colored as cyan and purple, respectively, (b) shows NNN paths that cannot be reached through one NN bonds. Yellow and grey 
arrows represent paths connecting same (Ir4 to Ir4) and different sublattices (Ir4 to Irl in the figure), respectively, (c) shows all third NN 
hopping paths. Purple and red arrows show bonds connecting different and same sublattices, respectively. Note that, bonds in (b) and red 
bonds in (c) does not have any counterparts in the 2D honeycomb lattice. 


t 2 g hopping terms and tight-binding bands from Wannier 
orbitals 

Table S3 shows a partial list of Ir t 2 g hopping terms (up to 
third NN) from the Wannier orbitals, where the convention for 
the coordinate system and the illustration of NNN and third 
NN hopping terms are in Fig. S2. Full list of hopping terms 
can be restored by applying the C 2 rotations and inversion 
operations at the centers of Z and X (X’) bonds, respectively. 
Three C 2 rotations — C 2 — are allowed, where a = x-\-y 
and b = y — X. 

Contrary to the relatively simple NN hopping channels as 
shown in the inset of Fig. S2(a), a number of distinct NNN 
hopping terms are introduced due to the three-dimensional 
twisting of the honeycomb lattice [42]. The NNN hopping 
channels can be classified into two kinds, depending on 


whether they are analogous to the NNN hopping in the 2D 
honeycomb lattice or not. Fig. S2(b) shows the 2D-like NNN 
hoping channels, which can be reached through one interme¬ 
diate NN Ir site. Depending on whether they belong to same 
zigzag chain composed of only X (or X’) bonds or connect 
different chains, they are divided into two different classes 
^NNN- Hopping amplitudes in these channels are 
larger than the other channels, and which can¬ 

not be reached through one Ir site as shown in Fig. S2(c). 
There are also non-negligible third NN hopping terms, 
and t 3 NN’ which can be seen in Fig. S2(d) and Table S3. Like 
NNN hopping channels, third NN channels can be classified 
depending on whether they have their 2D counterparts or not. 
4 nn resembles the third NN hopping channel in the 2D hon¬ 
eycomb lattice, while is similar to the interlayer hopping 
terms in aAlO series. 
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FIG. S3. (Color online) Band structures (solid red lines) calculated from tight-binding calculations including up to (a) NN, (b) NNN, and (c) 
third-NN hopping terms from the Wannier orbital calcualtions tabulated in Table. S3, compared to those from ab-initio results (dashed grey 
lines). SOC is included in the calculations. 


The role of further-neighbor hopping amplitudes in the 
band structure are shown in Fig. S3, where the evolution of 
band structure as we include NNN and third NN terms are 
presented. Fig. S3 shows the change of the band structures 
as further-neighbor hopping terms are included with the pres¬ 
ence of SOC; Fig. S3 (a), (b) and (c) show the bands with hop¬ 
ping terms up to NN, NNN, and third NN terms, respectively, 
with the presence of SOC. One can see that, the large SOC 
in Ir tends to make the jeff = 1/2 bands to be flatter in this 
locally honeycomb-like lattice, and including NNN and third 
NN terms does not change the overall behavior. Comparing 
Fig. S3 (a) and Fig. S3(b), one can notice that the dispersion 
inside the jeff =1/2 subbands is affected by the NNN terms, 
but the semi-metallic character is left unchanged. Inclusion 
of third NN terms, as can be seen in Fig. S3(c), makes the 
dispersion slightly closer to the ab-initio bands. 


Supplementary Material D: 


NN exchange interactions 


The exchanges J, K, and F are given by (suppressing the 
bond label a) 


4 

21 [ m (2f/ + 3A)2 


K = 


^2Jh r {ti — —3^2 

9 [ (2f/ + 3A)2 


,r = 


64Ji^ t2{ti — ts) 

9 (2f/ + 3A)2’ 


( 1 ) 


where U {i = 1,2,3), Jh, U, and A are the NN hopping 
amplitudes, Hund’s coupling, on-site Coulomb repulsion, and 
SOC respectively [9]. U is illustrated in Fig. S2(a). Note that, 
the small amount of NN Heisenberg interaction is attributed 
to the cancelation between the 2ti and ts in the antiferromag¬ 
netic contrubution to J in Eq. I. Since t 2 is the largest term, 
as mentioned in the main text, ferromagnetic K becomes the 
most dominant contribution in the exchange interactions. 


TABLE S3: A subset of Ir t 2 g hopping terms Tij as representatives of each 
hopping channels up to third NN, where Hhop = 

being the creation and annihilation operator for t 2 g states, respectively, d is 
approximate distance between Ir and O. Other hopping terms can be recovered 
by applying Tji = T^, rotations, and inversion operations. 


Kind 

Tij (in Cartesian coord.) 

Sublattice 



t/eff = 

: 0.0 eV 


f7eff = 

: 1.5 eV 


t/eff 

= 3.0 eV 

tNN 





dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

X,X’ 

(-d, 0,-\-d) 

1 - 

4 

dxy 

-hO.088 

-hO.018 

+0260 

+0.0S0 

-hO.019 

+0.216 

-hO.064 

+0.021 

+0.2S9 





dxz 

-hO.018 

-0.152 

-hO.013 

-hO.020 

-0.110 

+0.013 

-hO.021 

-0.051 

0.005 





dyz 

-hO.259 

-hO.013 

-hO.078 

+0216 

+0.013 

+0.061 

-hO.288 

0.003 

+0.052 

tNN 





dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

z 

(+d,-\-d, 0) 

1 - 

2 

dxy 

-0.162 

-0.022 

-hO.021 

-0.119 

-0.024 

+0.023 

-0.059 

-0.031 

-hO.030 





dxz 

-hO.016 

-hO.087 

-0.239 

+0.011 

-hO.078 

-0.255 

-hO.025 

+0.012 

-0.269 





dyz 

-0.016 

-0.239 

-hO.086 

-0.017 

-0.254 

+0.011 

-0.024 

-0.271 

+0.056 

Continued in next page... 
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Continued from previous page... 


Kind 

Vij (in Cartesian coord.) 

Sublattice 



t/eff = 

: 0.0 eV 


Vff = 

: 1.5 eV 


t/eff 

= 3.0 eV 

^NNN 




dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

dxy 

dxz 

dyz 


(- i -( i ,+ 2 d ,-( i ) 

1 ^ 3 

dxy 

0.002 

- 0.012 

+0.039 

0.001 

-0.015 

-hO.044 

0.001 

-0.018 

+0.041 




dxz 

+0.013 

0.001 

-hO.Oll 

-hO.018 

0.001 

+0.014 

+0.024 

0.001 

+0.011 




dyz 

+0.063 

0.004 

0.002 

+0.015 

0.007 

0.000 

-hO.089 

- 0.010 

0.001 

%NN 




dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

dxy 

dxz 

dyz 


(- d ,- i -( i ,- i - 2 ( i ) 

1 ^ 1 

dxy 

0.002 

0.008 

-0.014 

0.003 

- 0.011 

-0.017 

0.003 

-0.014 

- 0.020 




dxz 

+0.014 

0.001 

+0.039 

+0.011 

0.000 

+0.045 

+0.020 

0.002 

-hO.050 




dyz 

0.008 

+0.015 

0.001 

+0.011 

-hO.089 

0.000 

+0.014 

-hO.103 

0.001 

/III 

%NN 




dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

dxy 

dxz 

dyz 


(- d ,- h ( i ,- 2 ti ) 

1 ^ 1 

dxy 

0.001 

-hO.Oll 

0.001 

0.000 

-hO.013 

0.000 

0.002 

-hO.015 

0.003 




dxz 

0.001 

0.007 

-hO.038 

0.000 

0.005 

+0.045 

0.003 

0.002 

+0.051 




dyz 

- 0.011 

-hO.030 

0.008 

-0.013 

+0.036 

0.006 

-0.015 

+0.041 

0.004 

/IV 

%NN 




dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

dxy 

dxz 

dyz 


{+d,-2id,+d) 

1^4 

dxy 

-hO.012 

0.007 

-0.030 

0.009 

0.008 

-0.035 

0.005 

0.008 

-0.041 




dxz 

0.007 

0.008 

0.006 

0.008 

0.009 

0.008 

0.008 

+0.011 

0.008 




dyz 

-0.030 

0.007 

0.003 

-0.035 

0.008 

0.007 

-0.041 

0.008 

- 0.012 

4 nn 




dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

dxy 

dxz 

dyz 


( 0,+2dr2d) 

1 ^ 2 

dxy 

0.007 

-0.013 

-0.014 

0.001 

-0.014 

-0.015 

0.008 

-0.015 

-0.017 




dxz 

0.008 

0.006 

0.006 

0.009 

0.008 

0.007 

- 0.011 

+0.011 

0.009 




dyz 

-0.014 

-0.016 

-0.034 

-0.015 

-0.018 

-0.030 

-0.017 

- 0.020 

- 0.022 

/II 

^3NN 




dxy 

dxz 

dyz 

dxy 

dxz 

dyz 

dxy 

dxz 

dyz 


(-d,+d, 0 ) 

1 ^ 1 

dxy 

-0.045 

0.006 

- 0.012 

-0.046 

0.001 

-0.013 

-0.046 

- 0.010 

-0.013 




dxz 

+0.012 

0.008 

0.001 

+0.013 

0.009 

0.001 

+0.013 

+0.012 

0.002 




dyz 

0.006 

-0.015 

0.008 

0.008 

-0.017 

0.009 

+0.010 

-0.018 

+0.012 
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